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ABSTRACT 


Model experiments on Alfvén’s theory of aurorae, performed earlier by the author, resulted in 
auroral zones with a smaller eccentricity than theoretically expected if space charges are neglected. 
Preliminary measurements of the charges in the model have now been made, and on the basis of 
this it is possible to theoretically account for the smaller eccentricity. From the theory it also 
comes out that the measured space charge distribution will tend to increase the polar distance of 
the auroral zones. The influence of variable terrella voltage has also been investigated and the 
result is a good agreement between theory and experiment. 


Introduction 


Model experiments on Alfvén’s theory of aurorae (Alfvén, 1950) have been per- 
formed earlier by Block (1955). These resulted in a lower eccentricity than indi- 
cated by the theory, which neglects space charges. It has been suggested that space 
charges may account for the difference (Block, 1956). Another suggestion was that 
the difference may partly be accounted for by the fact that the terrella potential in 
the experiments differed considerably from that assumed in the theory (Block, 1955). 
This paper presents some results of measurements of space charges and of the in- 
fluence of the terrella potential together with theoretical calculations which largely 
agree with the experimental results. 


Experiments 


Measurements of the space charges in the model have been attempted by G. 
Apelqvist and C. G. Fialthammar (unpublished). They put probes in the discharge 
region (dark discharge) and measured the potential distribution. Their measurements 
could not be fulfilled because of an accident with the apparatus, so we do not know 
the potential in very much detail. However, sufficient material is available to get a 
preliminary rough picture of the charge distribution in the equatorial plane. Since 
this agrees with what we should qualitatively expect we publish it here with the 
reservation that more complete measurements will modify it to some extent. 

The result is shown in Fig 1. The upper part of the figure shows the difference, 
AV, between the potential along the x-axis in the model when there is a discharge 
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Fig. 1. Space charge distribution 
in the equatorial plane of the 
experimental discharge and its in- 
fluence on the potential along the 
@-axis, which is directed in the 6- 
hour direction through the earth’s 
centre. V,/V,=0.5 which gives no 
net current to the terrella. 


and the potential when there is no discharge but unchanged potentials on all elec- 
trodes, i.e. the potential due to space charges as a function of x. It is seen that AV is 
always positive but less positive on the cathode side of the terrella than on the anode 
side. The lower part of the figure shows the sign of the space charges in the equatorial 
plane. The distribution of charges differs from that in Alfvén’s auroral theory, 
because it is turned 90° to the “‘east’’ so that we have negative space charge on the 
cathode side. However, this agrees with what should be expected, because as pointed 
out by Block (1956) the magnetic field is too weak to make the radius of curvature 
of the ion orbits small in comparison to the forbidden zone, so the ions move essentially 
along the electric field. Hence the anode side of the terrella acts as a collector of 
ions and the cathode side as an anode, collecting electrons. The total current to the 
terrella is zero, but there is a current from the anode side to the cathode side (also 
observed in the experiments), corresponding to a current from the evening side to 
the morning side in the auroral zones of the earth as explained in an earlier paper 
(Block, 1956, p. 317). According to the theory of probes in discharges, there should 
be positive space charge outside a probe collecting ions, and negative space charge 
outside a probe collecting electrons. Thus, if the terrella is regarded as a sort of 
double probe the above picture agrees with probe theory. 

The influence of the terrella potential on the location of the auroral zones has 
been investigated in the dark discharge. The potential difference, Vz, between anode 
and cathode was kept constant, V;=5 kV, the magnetic field at the poles of the 
terrella B, = 1000 gauss, the pressure p = 8 x 10-4 mm Hg and the discharge current 
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I, =0.6 mA. The potential, V,, of the anode relative to the terrella was varied. The 
result of the measurements is shown in Fig. 2. 6, is the maximum angular polar 
distance and 6, is the minimum polar distance of the auroral zones. If V,/V_z=0 
the terrella is at anode potential, and if V,/V_z=1 it is at cathode potential. 

It is seen from Fig. 2 that the more V/V, increases, the closer do the auroral zones 
come to the poles and the larger is the eccentricity. 


Theory of the influence of the terrella potential 


The equations of the electron orbits in a magnetic dipole field with the dipole in 
the —z-direction in the presence also of a homogeneous electric field Hy in the x- 
direction and a homogeneous magnetic field By in the z-direction has been obtained 
by Alfvén (loc. cit.). If the terrella is charged, however, the electric field is not homo- 
geneous any more, and then Alfvén’s equations are modified in the following way 
where the symbols have their usual meaning (cf. Alfvén, 1950, p. 179, and Block, 
1955, p. 85). 


V=V,+ Ey (tp — 2) + O/R, (1) 
B=B,+/R', (2) 
pp =eV/B. (3) 


Q is a constant proportional to the charge of the terrella. The terrella has been 
assumed to be infinitely small so that the charge can be regarded as a point charge. 
However, in order to make comparisons with the experiments we must assume a 
finite terrella radius Ry, in which case the terrella potential in relation to the mean of 
the anode and cathode potential is 


V,=Q/R. (4) 

Combining (1), (2) and (3) we find 
Lp —x =I[A/R* — K/R, (5) 
IA = b/e E,; K =Q/E, = V,Ro/ Eo. (6) 
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Fig. 3. The forbidden zone curve in the 
equatorial plane. 


4 
If Q=0 we have tp=$V3L=1.76L, 
4 
a, =V3L=1.32 L (7) 
x, = —0.74L. 


For the definition of x, and 2, see Fig. 3. If @+ 0 we find | 


2,=V—-4K+V31A+4K2 ees 
with Up =x, + LA/x3 — K/ax, (9) 
and x, is determined by tp — ay = Kye Lx, (10) 


In Fig. 4 the results of (8) and (10) have been plotted as a function of V/V, where 
according to (4) and (6) V,/V, is determined by 


V,/Va=0.5 — V,/Vq=0.5 — K Ey/ Ro Va. (11) 


For reasons which will become clear in the next paragraph we have chosen L = 2.20 
Ry when calculating the theoretical curves in Fig. 4. (cf. eq. (27)). As a comparison 
the experimental values of x, and 2x,, obtained by projection of the experimental 
auroral zone along the dipole lines, are shown in the figure. It is seen that theory 
and experiment agree, in as much as that they both indicate increasing x, and a, with 
increasing V,/V4,, i.e. decreasing terrella potential. The rate of increase of 2, and 
X_ respectively are also about the same theoretically and experimentally. 

If we had chosen ZL < 2.20 R, the theoretical curves would have been lowered 
correspondingly but as will be seen in the next praragraph it is possible to account 


for the discrepancy between theory and experiment by considering the influence of 
space charges. 
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Fig. 4. The distances x, and «, as functions 
of the terrella voltage according to the theory 
neglecting space charges. The rings and crosses 
indicate experimental values. 


Theory of the influence of space charges 


We assume that the space charges change the potential by 


AV =AV(z,y,2), (12) 
where AV 0 at infinity. This means that the equations of the electron orbits are 
P= Ver hit, 2) A Pe, (13) 
together with (2) and (3). A combination of these gives 
tp —«2 = IA/R§—AV/Ey. (14) 
Differentiation with respect to 2 gives 
chante ge ea er ry (15) 


where R? = x? + y?. 
Putting y =0, yy’ =0 and (¢AV/éy)y’ =0, we obtain 
3 [4 1 dAV 


=—+— —. 16 
1 Ew az (16) 
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The solution of this equation is 2 = 29, 


{ 3 i 
cerry E,) (@AV/0x)) cnx, 
/ 


If AV=0 this goes over in the usual value of 
4 
Le =o = V3L. 
The corresponding 2p is found by putting «=a, y =90 in (14) 
Gp = Xy + L*/x8 — AV (xq)/ Eo 


where A V (x,) is the value of AV at x =a, y =z =0. Since 


a ier 


we get approximately — 1 
where Lp, =% V3.L 


is the usual value of x, for AV =0. 
Finally we find x, from (14) by putting R = —a,, x=2: 


Ly = Xp + L*/2x# A V (a) / Eo: 


; TA 
or using (19) Dp i me a {AV (x) — AV (a5)}. 
Uy Ey 


0 


i 1 
Introducing D= i {AV (x,) — AV (a)} 
0 
we find from (20a) 
@D/éx,=1+ oL/ a, 
which gives an approximate solution of (20a) 


D 


BL BE WEE 5 


Now the usual value of x, for AV=0 is 


10 i 0.74 fh 
Introducing this in (23) we get 
D 
% = Lg t il =Mot 7a oe {AV (x) — AV (x,)}. 
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According to the experimental curve in Fig. 1 AV is always positive but AV (z,) is 
slightly larger than A V (x,) so that 2, is slightly larger than x. This means that the 
auroral zone is pushed somewhat towards the equator here. At x = 2, it is seen from 
(17) that the same happens if dA V/éx <0. This appears to be the case in the model. 

Qualitatively the results of this treatment is that space charges seem to increase 
the polar distance of the auroral zone. This means that Z can be somewhat larger 
and the beam fields Hy) and By accordingly smaller than assumed earlier. 


Application to the experiments 


A numerical application of these results to the experiments is very uncertain be- 
cause of the still insufficient potential data available. However, an attempt will be 
made in order to see what is needed to explain the well established experimental data 
given in Fig. 4, The only potentials we need then is A V (x,) — AV (a,). This seems to 
be of the order of 100 volts. 

The total discharge voltage was 5 kV and the distance between cathode and 
anode is 11 Ry. Hence 


Ey, = 5000/11 Ry. (25) 
Then (24) and Fig. 4 (V,/V_=0.5) gives 
R 

a= —0.74L 2-100 = — 1.61 Ry. 26 

Ts 0.74 D+ 5000 ° (26) 

From this we get L = 2.20 Ro. (27) 


Fig. 4 shows that the experimental x, = 2.0 Rp. 
Equation (17) gives 
(@AV /0x),-2,= Ey (1—3 LA/23). (28) 


Our numerical values give 


4 5000 (: AS 1550 


f 2 ee 29 
Oa i1R, 4 R, ee) 


This means a field strength due to the space charges of about 300 V/cm in addition 
to #, which is only 85 V/em. In view of the curves shown in Fig. 1 this does not 
seem unreasonable. 

Measurements of the total current to the terrella revealed that this was zero when 
V./V_=0.5. Thus if V, is more negative the terrella mainly receives ions and the 
electrons are reflected away. It was also noticed that if V,/Vq > 0.55 a glow discharge 
was initiated. This may be due to intense ionization by oscillating electrons, reflected 
by the negative terrella potential. In any case the space charge should be radically 
changed when we approach this state. This may account for the sharp bend of the 
experimental curve for 2, shown in Fig. 4 and, which is the same, the strong increase 
of the eccentricity at V/V > 0.5 in Fig. 2. However, it may also be that the bend of 
the curve is not quite as sharp as shown in the figure, because the experimental errors 
in ¢ may be as large as one tenth of the measured values. According to our theory 
the eccentricity should be very insensitive to the terrella voltage if space charges 
are neglected. 
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_ From the experimental curves in Fig. 4 we find for V, /Vq=9.5 
4 (a, +2) = 1.80 Ro. 
According to the simplified theory which neglects space charges we have 
4 (x, + aq) = 1,05 0. 


This means that the simple theory gives the same auroral mean latitudes as our 


experiments if 
L=1.75 Ro. 


However, when we accounted for space charges we found in eq. (27) that L = 2,20 Ro. 
Hence, the more complete theory can give the same auroral mean latitudes as the 
simple theory provided L is increased by a factor 


k = 2,20/1.75 = 1.26. 


Since L ~ E)°”> the beam fields H, and B, may be only 40 % of what has been assumed 
before, and still produce aurorae at the same latitudes, if these experimental results 
are applicable also to the cosmical phenomena. This is, however, uncertain because 
the space charge distribution in nature is different from that in the model due to 
the relative difference in the radius of curvature of the ion orbits. On the other hand, 
it may be that this does not matter so much if only the electric field at x, is strongly 
increased by space charges in nature. This would increase the drift velocity due to 
the electric field so that the point x, where the magnetic gradient drift balances the 
electric field drift, is pushed closer to the earth. 

These results indicate very clearly that Alfvén’s auroral theory is applicable to 
the model experiments. 
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